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. How to measure the effective action for disordered systems 
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In contrast to standard critical phenomena, disordered systems need to be treated via tlie 
Functional Renormalization Group. Tlie latter leads to a coarse grained disorder landscape, 
which after a finite renormalization becomes non-analytic, thus overcoming the predictions 
of the seemingly exact dimensional reduction. We review recent progress on how the non- 
analytic effective action can be measured both in simulations and experiments, and confront 
, theory with numerical work. 
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1. Introduction 



When talking to his experimental colleagues about the marvels of field theory and 
\^ ' his recent achievements in computing the effective action of his favorite model, 

OO . the conversation is likely to resemble this: 

Theorist : I have a wonderful field theory, I can even calculate the effective action! 
Experimentalist: Can I see it in an experiment? Can I measure it? 



Theorist : . . . well, that's difficult, but it tells you all you want to know . . . 
' Experimentalist: okay, I understand, another of these unverifiable predictions, . . . 

■ Here we will see how to measure it, considering the explicit and far-from- 

trivial example of elastic manifolds in a disordered environment. Due to a lack of 
space, we will not be able to give all arguments in the necessary details. We rec- 
ommend that the reader consults the recent "Basic Recipes and Gourmet Dishes", ' 



I to which we also refer for a more complete list of references. 



2. The disordered systems treated here - our model 

Let us first give some physical realizations. The simplest one is an Ising magnet. 
Imposing boundary conditions with all spins up at the upper and all spins down 
at the lower boundary (see figure 1), at low temperature T, a domain wall will 
form in between. In a pure system at T = 0, this domain wall is completely 
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flat; it will be roughened by disorder. Two types of disorder are common: random 
bond (which on a course-grained level represents missing spins) and random field 
(coupling of the spins to an external random magnetic field). Figure 1 shows, how 
the domain wall is described by a displacement field u{x). Another example is the 
contact line of water (or liquid hydrogen), wetting a rough substrate. A realization 
with a 2-parameter field u{x) is the deformation of a vortex lattice; the position of 
each vortex is deformed from xto x + u{x). A 3-dimensional example are charge 
density waves. 

All these models are described by a displacement field 

xeM.'^ — > u{x) e . (1) 

For simplicity, we now set iV = 1. After some initial coarse-graining, the energy 
W = T-Lc\ + Wdo consists out of two parts: the elastic energy, and the disorder: 

HeiM = j d'^xi (Vw(x))' , HdoM = j d''xV{x,u{x)) . (2) 

In order to proceed, we need to specify the correlations of disorder:-' 

V{x,u)V{x',u') 5'^{x - x')R{u - u) . (3) 
Fluctuations u in the transversal direction will scale as 

[u{x) - u{y)f ^ \x - yf'^ . (4) 

There are several useful observables. We already introduced the roughness- 
exponent C,. The second is the renormalized (effective) disorder function R{u), 
and it is this object we want to measure here. Introducing replicas and averaging 
over disorder, we can write down the bare action or replica-Hamiltonian 

^M = ^E / d'a;i(Vu,(x))'-^ / d'' x R{ua{x) - u,{x)) . (5) 



u(x) 





defect 

("random bond") 



' 'random field' ' 




Fig. 1. An Ising magnet at low temperatures forms a domain wall described by a function u{x) 
(right). An experiment on a thin Cobalt film (left);- with kind permission of the authors. 
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Let us stress that one could alternatively pursue a dynamic or a supersymmetric 
formulation. Since our treatment is perturbative in R{u), the result is unchanged. 

3. Dimensional reduction 

There is a beautiful and rather mind-boggling theorem relating disordered systems 
to pure systems (i.e. without disorder), which applies to a large class of systems, 
e.g. random-field systems and elastic manifolds in disorder. It is called dimen- 
sional reduction and reads as follows:"^ 

Theorem: A d-dimensional disordered system at zero temperature is equivalent to 
all orders in perturbation theory to a pure system in d ^ 2 dimensions at finite 
temperature. 

Experimentally, one finds that this result is wrong, the question being why? Let 
us stress that there are no missing diagrams or any such thing, but that the problem 
is more fundamental: As we will see later, the proof makes assumptions, which 
are not satisfied. Before we try to understand why this is so and how to overcome 
it, let us give one more example. We know that the width u of a d-dimensional 
manifold at finite temperature in the absence of disorder scales as m x^^"'^^^'^. 
Making the dimensional shift implied by dimensional reduction leads to 

[u{x) - u{0)f ~ x^''^ = a;2« i.e. C = • (6) 

4. The Larkin-length 

To understand the failure of dimensional reduction, let us turn to an interesting 
argument given by Larkin.^ He considers a piece of an elastic manifold of size 
L. If the disorder has correlation length r, and characteristic potential energy /, 

d 

this piece will typically see a potential energy of strength -Edo = / (7)^ ■ On 
the other hand, there is an elastic energy, which scales like Ed ~ cL'^~^ . These 

jjr 1 . More 

important than this value is the observation that in all physically interesting dimen- 
sions d < 4, and at scales L > Lc, the membrane is pinned by disorder; whereas 
on small scales the elastic energy dominates. Since the disorder has a lot of min- 
ima which are far apart in configurational space but close in energy (metastability), 
the manifold can be in either of these minimas, and the ground-state is no longer 
unique. However exactly this is assumed in the proof of dimensional reduction. 

5. The functional renormaUzation group (FRG) 

Let us now discuss a way out of the dilemma: Larkin's argument suggests that 
d = 4 is the upper critical dimension. So we would like to make an e = 4 — d 
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expansion. On the other hand, dimensional reduction tells us that the roughness is 
( = (see (6)). Even though this is systematically wrong below four dimen- 
sions, it tells us correctly that at the critical dimension d = 4, where disorder is 
marginally relevant, the field u is dimensionless. This means that having identi- 
fied any relevant or marginal perturbation (as the disorder), we can find another 
such perturbation by adding more powers of the field. We can thus not restrict 
ourselves to keeping solely the first moments of the disorder, but have to keep 
the whole disorder-distribution function R{u). Thus we need afunctional renor- 
malization group treatment (FRG). Functional renormalization is an old idea, and 
can e.g. be found in.^ For disordered systems, it was first proposed in 1986 by D. 
Fisher.^ Performing an infinitesimal renormalization, i.e. integrating over a mo- 
mentum shell a la Wilson, leads to the flow deR{u), with (e = 4 — d) 

deR{u) = (e - 4C) R{u) + CuR'{u) + ^R"{uf - R"{u)R"{Q) . (7) 

The first two terms come from the rescaling of R and u respectively. The last two 
terms are the result of the 1-loop calculations, see e.g.'-^ 

More important than the form of this equation is it actual solution, sketched in 
figure 2. After some finite renormalization, the second derivative of the disorder 
R"{u) acquires a cusp at u = 0; the length at which this happens is the Larkin- 
length. How does this overcome dimensional reduction? To understand this, it is 
interesting to study the flow of the second and forth moment. Taking derivatives 
of (7) w.r.t. u and setting u to 0, we obtain 

diR"{{)) = (e - 20 i?"(0) + R"'{Qf [e - 2C) R"{0) (8) 
diR""{Q) = ei?""(0) + 3i?""(0)2 + 4i?"'(0)i?""'(0) ei?""(0) + iR""{Qf. 

Since R{u) is an even function, and moreover the microscopic disorder is smooth 
(after some initial averaging, if necessary), R"'{Q) and R'""{0) are 0, which 
we have already indicated. The above equations for R"{0) and R""{0) are in 
fact closed. The first tells us first that the flow of R"{0) is trivial and that 
= e/2 = This is exactly the result predicted by dimensional reduction. 
The appearance of the cusp can be inferred from the second one. Its solution is 
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Fig. 3. Generation of the cusp, as explained in tlie main text. 

i?""(0)|^ = i_3e(e°'/_i)/, , with c = R""{0)l^Q- Thus after a finite renormaHza- 
tion i?""(0) becomes infinite: The cusp appears. By analyzing the solution of the 
flow-equation (7), one also finds that beyond the Larkin-length i?"(0) is no longer 
given by (8) with i?"'(0)^ = 0. The correct interpretation of (8), which remains 
valid after the cusp-formation, is deR"{Q) = (e - 2C) R"{0) + i?"'(0+)2. Renor- 
malization of the whole function thus overcomes dimensional reduction. The ap- 
pearance of the cusp also explains why dimensional reduction breaks down: The 
simplest way to see this is by redoing the proof for elastic manifolds in disorder, 
which in the absence of disorder is a simple Gaussian theory. Terms contributing 
to the 2-point function involve i?"(0), TR""{0) and higher derivatives of R{u) at 
u = 0, which all come with higher powers of T. To obtain the limit of T ^ 0, 
one sets T = 0, and only i?"(0) remains. This is the dimensional-reduction result. 
However we just saw that R""{Q) becomes infinite. Not surprisingly R""{0)T 
may also contributes; indeed one can show that it does, hence the proof fails. 

6. The cusp and shocks 

Let us give a simple argument of why a cusp is a physical necessity, and not an 
artifact. The argument is quite old and appeared probably first in the treatment 
of correlation-functions by shocks in Burgers turbulence. It was nicely illustrated 
in.^ Suppose, we want to integrate out a single degree of freedom coupled with 
a spring. This harmonic potential and the disorder term are represented by the 
parabola and the lowest curve on figure 3(a) respectively; their sum is the re- 
maining curve. For a given disorder realization, the minimum of the potential as a 
function of u is reported on figure 3(b). Note that it has non-analytic points, which 
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mark the transition from one minimum to another. Taking the derivative of the po- 
tential leads to the force in figure 3(c). It is characterized by almost linear pieces, 
and shocks (i.e. jumps). Calculating the force-force correlator, the dominant con- 
tribution for small distances is due to shocks. Their contribution is proportional to 
their probability, i.e. to the distance between the two observable points. This leads 
to F{u)F{Q) — F(0)2 — c\u\, with some numerical coefficient c. 

7. The field-theoretic version 

The above toy model can be generaUzed to the field theory.' Consider an interface 
in a random potential, and add a quadratic potential well, centered around w: 

r irn ^ 

wrotM = j —Hx) - wf + KiM + HdoN ■ (9) 

In each sample (i.e. disorder configuration), and given w, one finds the minimum 
energy configuration. This ground state energy is 

V{w):=Tmn'HZAn\- (10) 
It varies with w as well as from sample to sample. Its second cumulant 

V{w)V{w') = L'^R{w ~w') (11) 

defines a function R{w) which is proven** to be the same function computed in 
the field theory, defined from the zero-momentum effective action."^ Physically, 
the role of the well is to forbid the interface to wander off to infinity. The limit of 
small m is taken to reach the universal limit. The factor of volume is neces- 
sary, since the width of the interface in the well cannot grow much more than 
m^^. This means that the interface is made of roughly i/L,„ pieces of internal 
size L,n ~ m pinned independently: (11) expresses the central-limit theorem and 
R{'w) measures the second cumulant of the disorder seen by each piece. 

The nice thing about (11) is that it can be measured. One varies w and 
computes (numerically) the new ground-state energy; finallying averaging over 
many realizations. This has been performed recently in" using a powerful exact- 
minimization algorithm, which finds the ground state in a time polynomial in the 
system size. In fact, what was measured there are the fluctuations of the center of 
mass of the interface u{w) = L^'^ J d'^x uo{x; w): 

[w - u{w)][w' - u(w')]'' = m-^L-'^A{w - w') (12) 

which measures directly the correlator of the pinning force A(u) = —R"{u). 
To see why it is the total force, write the equilibrium condition for the center of 

mass w? [w — u[w)] + L^"^ J d'^x F{x, u) = (the elastic term vanishes if we use 
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z (z/4forRB) 

Fig. 4. Filled symbols show numerical results for Y{z), a normalized form of the interface displace- 
ment correlator -~R"{u) [Eq. (12)], for D = 2 + 1 random field (RF) and D = 3 + 1 random 
bond (RB) disorders. These suggest a hnear cusp. The inset plots the numerical derivative Y' {z), with 
intercept Y' (0) Ri —0.807 from a quadratic fit (dashed line). Open symbols plot the cross-coiTelator 
ratio Ys (z) = A12 (z)/ An (0) between two related copies of RF disorder. It does not exhibit a cusp. 
The points are for confining wells with width given by A/2 = 0.02. Comparisons to 1-loop FRG 
predictions (curves) are made with no adjustable parameters. Reprinted from. ' ' 

periodic b.c). The result is represented in figure 4. It is most convenient to plot the 
function Y = A(u)/A(0) and normalize the u-axis to eliminate all non-universal 
scales. The plot in figure 4 is free of any parameter. It has several remarkable 
features. First, it clearly shows that a linear cusp exists in any dimension. Next it 
is very close to the 1-loop prediction. Even more remarkably the statistics is good 
enough" to reliably compare the deviations to the 2-loop predictions of.'^ 

When we vary the position w of the center of the well, it is not a real motion. 
It means to find the new ground state for each w. Literally "moving" w is another 
interesting possibility: It measures the universal properties of the so-called "de- 
pinning transition".'^ ''* This was recently implemented numerically (see Fig. 6). 




Fig. 5. The measured Y{u) with the 1- and 2-loop coiTections subtracted. Left: RB-disorder, right: 
RF-disorder. One sees that the 2-loop corrections improve the precision. 
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Fig. 6. Running tlie RG in a numerical simulation: Crossover from RB disorder to RF for a driven 
particle (left)." Residual eiTor for Y{u) for a driven string'- which show that statics and depinning 
are controled by different fixed points. 
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